Transverse force on a vortex and vortex mass: 
effects of free bulk and vortex-core bound quasiparticles 
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The paper reassesses the old but still controversial problem of the transverse force on a vortex 
and the vortex mass. The transverse force from free bulk quasiparticles on the vortex, both in 
the Bose and the Fermi liquid, originates from the Aharonov-Bohm effect. However, in the Fermi 
liquid one should take into account peculiarities of the Aharonov-Bohm effect for BCS quasiparticles 
described by two-component spinor wave functions. There is no connection between the transverse 
force (either from free bulk quasiparticles or from vortex-core bound quasiparticles) and the spectral 
flow in the vortex core in superfluid Fermi liquid, in contrast to widely known claims. In fact, there 
is no steady spectral flow in the core of the moving vortex, and the analogy with the Andreev bound 
states in the SNS junction, where the spectral flow is really possible, is not valid in this respect. 

The role of the backflow on the vortex mass is clarified. The backflow is an inevitable consequence 
of a mismatch between the currents inside and outside the vortex core and restores the conserva- 
tion of the particle number (charge) violated by this mismatch. In the Fermi liquid the backflow 
compensates the current through the core bound states, which is a source of the vortex mass (the 
Kopnin mass). This results in renormalization of the Kopnin vortex mass by a numerical factor. 

PACS numbers: 67.25. dk,67.30.he,74.25.Uv 
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I. INTRODUCTION 

Discussions and debates on the transverse force on 
a vortex in superfluids (neutral and charged) continue 
during many decades and have been a topic of reviews 
and booksi~— . They focused on quasiparticle contribu- 
tions to this force, which are connected with geomet- 
rical (Aharonov-Bohm-Berry) phases in the superfluid 
around the vortex. Although in most of practical cases 
the vortex can be considered as a massless object gov- 
erned by the gyroscopic dynamics, the concept and the 
magnitude of the vortex mass was also vividly discussed 
in the Bose and the Fermi superfluids^ - — . 

Despite a huge literature on this subject, there still re- 
main some issues, which require further clarification, es- 
pecially for the Fermi superfluids. The present paper ad- 
dresses these issues. One of them is the concept of spec- 
tral flow, which is sometimes assumed^^ to be a mech- 
anism responsible for the transverse force produced by 
bound states in the vortex core (Kopnin-Kravtsov force) 
and for the part of the transverse force from scattering 
of free quasiparticles in the Fermi superfluid. This paper 
argues that the spectral flow cannot be responsible for 
any part of the transverse force simply because it is ab- 
sent in a core of a moving vortex, as already was noticed 
by Stone— in the past. On the other hand, all kinds of 
transverse forces can be understood within common ap- 
proaches like the scattering theory and the partial-wave 
expansion without any reference to the spectral flow. In 
particular, the part of the transverse force from scatter- 
ing of free bulk quasiparticles in the Fermi superfluid, 
which was presumed to originate from the spectral flow, 
directly follows from peculiarities of the Aharonov-Bohm 
effect for BCS quasiparticles described by two-component 
spinor wave functions. 



Addressing the vortex mass, the present paper revises 
different contributions to it, compares them, and demon- 
strates what are their possible effects on vortex motion. 
In particular, it demonstrates that the so-called backflow 
vortex mass accompanies any type of vortex mass. The 
backflow mass is related with the kinetic energy of a su- 
perflow around the vortex core, which inevitably appears 
any time when the current density inside the core differs 
from that outside the core, and the intensity of the back- 
flow is determined from the continuity of the total fluid 
current. 

Let us present the nomenclature of various forces, 
which enter the equation of vortex motion: 



mn s K[z x (v L - v s ) = F n + F c 



The left-hand side is the Magnus force, which transfers 
momentum between the superfluid and the vortex. Here 
Vl and v s are the vortex and the superfluid velocities, 
m is the particle mass, n s is the superfluid density, and 
K = h/m is the circulation quantum. The force 



F n = -D(v L - v n ) - D'[z x (v L - v n )} 



(2) 



transfers momentum between the normal component (the 
gas of free quasiparticles) and the vortex. Here v n is the 
normal velocity. The coefficients D and D', 
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are determined by the longitudinal (transport) and the 
transverse cross-sections <7ii and a±, which will be de- 
termined further in the paper. Here /o(e) is the equi- 
librium Fermi distribution function of energy e of free 
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quasiparticles, and vq is the projection of the quasiparti- 
cle group velocity on the plane normal to the vortex line. 
The transverse force proportional to D 1 is the Iordanskii 
force. 

The force F c transfers momentum from the quasipar- 
ticles occupying bound states in the vortex core to im- 
purities in superconductors or to free bulk quasiparticles 
constituting the bulk normal component of the 3 He su- 
pcrfluid. The force has also two components, longitudinal 
and transverse to the relative normal velocity v n — vl, 
the latter called the Kopnin-Kravtsov forced. 

Finally, dP/dt is the inertial force, which is a product 
of the vortex mass and the vortex acceleration dv^/dt, 
the momentum P being the momentum of the vortex 
dependent on v^. 

The theory presented in this paper assumes that the 
quasiparticle mean-free path is much longer that the core 
size and therefore it cannot be used for high tempera- 
tures. The whole paper addresses neutral superfluids, 
although the results for the Fermi superfluids are rel- 
evant also for type II s-wave superconductors, since the 
effects of magnetic fields usually are not essential for vor- 
tex dynamics 2 . In superconductors the normal velocity 
v n usually vanishes in the coordinate frame related to 
the crystal lattice. 

The paper starts from Sec. [H]reminding the old results 
for semiclassical scattering of quasiparticles by a vortex. 
This shows the connection of the transverse force with 
the Aharonov-Bohm effect for quasiparticles. Section ITO1 
considers the scattering of BCS quasiparticles on the ba- 
sis of the Bogolyubov-de Gennes equations. The anal- 
ysis is done using the geometric optics and the partial- 
wave method. It demonstrates that the whole transverse 
force from free bulk quasiparticles is fully explained by 
the Aharonov-Bohm effect without referring to the con- 
cept of spectral flow. But one must take into account 
the peculiarities of the Aharonov-Bohm effect for BCS 
quasiparticles described by two- component spinor wave 
functions. Section HVl reminds properties of bound states 
in the vortex core in the Fermi superfluid focusing on the 
role of superfluid motion outside the core. Sections IVland 
IVII consider various contributions to the vortex mass in 
the Bose and the Fermi liquid respectively. Section IVIII 
discusses the derivation of the transverse force and the 
vortex mass from the Boltzmann equation focusing on 
the effect of superfluid transport past the vortex and on 
the comparison with the analysis of the previous sections. 
Section I Villi analyses possible effects of the vortex mass 
on vortex dynamics. Concluding discussion of the results 
is presented in Sec. IIXI Two appendices address more 
special issues: the simplified derivation of the spectrum 
of bound states for a core with linear growth of the gap 
as a function of the distance from the axis (App. [SJ and 
the derivation of the vortex mass for a core with linear 
growth of density in the Bose superfluid (App. |B|) . 



II. TRANSVERSE FORCE FROM THE 
SEMICLASSICAL SCATTERING THEORY 
(GEOMETRIC OPTICS) 

It is useful to start from the simplest approach to 
this problem based on the semiclassical scattering the- 
ory, which was first used for rotons by Lifshitz and 
Pitaevskii— long ago. 

The theory is based on the geometric optics. A quasi- 
particle moves along a well-defined trajectory and its mo- 
tion is described by variation of the position vector R 
and the momentum p of the quasiparticle in time. The 
classical Hamilton equations for them are: 



dR 
~~dt 



de 
dp' 



dp 

~dl 



dR' 



Here 



e(p) = e o(p) +P ■ v v 



(4) 



(5) 



is the energy of the quasiparticle in the moving fluid, eo 
is the quasiparticle energy in the resting fluid, and v v is 
the velocity induced by a rectilinear vortex: 



(6) 



2irr 2 



where r is a position vector in the plane normal to the 
vortex line (the projection of R on that plane). In order 
to simplify discussion we assume that the quasiparticle 
moves in the normal plane, so its momentum p lies in 
this plane. 

The vortex velocity field produces a force V(p • v v ) 
on the quasiparticle. The force may be considered as 
weak and the quasiparticle trajectory as nearly rectilin- 
ear. Suppose that the trajectory is parallel to the y-axis 
(Fig.fT^) and its impact parameter (the distance between 
the vortex line and the trajectory) is b — x. Then Eq. (|4]) 



gives 



dy 
dt 



~T7 = VG, 



V v ). 



(7) 



Here vq = de (p)/dp is the quasiparticle group veloc- 
ity in the resting fluid, which is in our case approxi- 
mately parallel to the axis y. Excluding time from these 
equations one has a differential equation determining the 
quasiparticle momentum variation along the trajectory: 



T = V(p ■«*,). 

dy v G 



(8) 



Integration of this equation assuming that the group ve- 
locity vg does not vary along the trajectory yields 



P(y) =P~ —v v (b,y), 

VG 



(9) 



where p = p(—oo) is the momentum at y — — oo. 

The scattering angle <p between the final and the initial 
momenta of the quasiparticle determines the momenta 
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FIG. 1. (Color online) Semiclassical scattering of rotons. a) 
Rotons at trajectories with x < or x > b* move past the vor- 
tex. Rotons at trajectories with < x < b* are fully reflected 
by the vortex. The shaded area (Andreev shadow) is classi- 
cally forbidden for rotons. b) Action variation 8S(b) along 
the trajectory as a function of the impact parameter b (di- 
mensionless variables), c) Scattering angle <p(b) as a function 
of the impact parameter b (dimensionless variables). 



p(l — cos ip) and p sin ip, which are longitudinal and trans- 
verse with the respect to the incident momentum p. The 
momenta are transferred by the scattered quasiparticle 
to the vortex. Correspondingly, the longitudinal and 
the transverse forces on the vortex from quasiparticles 
[see Eqs. @ and ©] are determined by the longitudinal 
(transport), 



a(ip)(l — cos (p)dip 



<P(bf 



db, 



(10) 



and the transverse, 

a± = a (ip) sin tp dip fa I ip(b)db, (11) 



effective cross-sections. Here 



cr(p) = 



db 

dp 



(12) 



is the differential cross-section. It was assumed that the 
scattering angle p « Px/p is small. 

In the Hamilton-Jacobi theory the momentum is con- 
nected with the classical action: p = dS/dr. Then 
Px — dSS(b)/db, where 



f°° p f°° 
SS(b) = / \p(y)-p]dy = / V y v vy 

J — oo J —oo 



dy. (13) 



is the variation of the classical action along the trajec- 
tory, which is a function of the impact parameter b. This 
yields: 



dSS(b) 
db 



db 



SS(-oo) - 5S(+oo) 



P 



(14) 



Bearing in mind that the velocity induced by the vortex 
is v v = «;V0(r) where the phase (j> = arctan (y/a;) is the 
azimuthal angle for the two-dimensional position vector 
r [see Eq. (|6])] one obtains that 



5S(b) 



pn 
2ttvq 



b 2 



-dy = — sign6 



pn 
2«g' 



(15) 



Eventually Eq. P^|) yields the transverse cross-section^ 



Vg 



(16) 



So we have obtained for the transverse cross-section an 
universal expression, which looks valid for any quasipar- 
ticle spectrum. The cross-section is proportional to the 
total variation of the classical action around the vortex 
line. Because of correspondence of the classical action to 
the quantum mechanical phase, this points out connec- 
tion of the transverse force with the geometric phase, or 
the Aharonov-Bohm effect. Equation (fT6|) yields a cor- 
rect transverse cross-section for phonons even though the 
semiclassical theory is not valid for phonons: there is no 
well-defined classical trajectories for phonons except for 
large impact parameters b at which the scattering angle 
tp is negligible. In the case of phonons the group velocity 
vg is the sound velocity c s . 

The simple universal expression (|14[) for the transverse 
cross-section, which does not depend on details of depen- 
dence of the scattering angle on the impact parameter, 
was obtained because an integrand in Eq. (fl"4|) is a deriva- 
tive of the action. But it does assume that the action is a 
continuous differentiable function of the impact parame- 
ter. Now we shall check it for rotons in superfiuid 4 He. 

The energy spectrum for rotons is eo(p) = A + (p — 
Po) 2 /2p, where A is the roton gap and p is the roton 
mass. According to the energy conservation law following 
from the Hamilton equations, Eq. (U), one has: 



A + 



[p{y) - pq] 

2/i 



■p{y)vv{y) 



(p - Po) 2 

2pL 



(17) 
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The right-hand side is the energy far from the vortex 
line, where p = p(— oo). The variation of the roton group 
velocity along the trajectory with the impact parameter 
b is given by 



vg(v) = p(yV \ Po = - x l (p - p ) 2 - 2pv vy 



where 




bb* 



b 2 + y 2 



(18) 



Here the characteristic scattering length 



npp 



^{p-Pof 



(19) 



is introduced and the vq = vg{— oo) = (p — Po)/ p is 
the roton group velocity far from the vortex line. The 
asymptotic expression Eq. (|9]) obtained at constant vq is 
valid for large impact parameters |6| 3> b* . 

In the classical scattering theory the point y — on the 
trajectory is a turning point: At y < the quasiparticle 
approaches to the scattering center (vortex line in our 
case) while at y > the quasiparticles moves away from 
the vortex line. Equation (|18[) shows that for impact 
parameters < b < b* the quasiparticle cannot reach 
the usual turning point y = since at y = —y*, where 
y* = y/b*b — b 2 , the group velocity changes a sign, and 
the quasiparticle starts to move back to y = — oo without 
an essential change of its momentum. This is Andreev 
reflection well known in the theory of superconductivity. 
At the point y = —y* p = po and the transition be- 
tween two branches of the roton spectrum with p > p 
(positive branch, parallel momentum and group velocity) 
and p < po (negative branch, antiparallel momentum and 
group velocity) occurs. Due to the Andreev reflection the 
shadow region is formed near the vortex line which is not 
available for the roton classical trajectories. That shadow 
(Andreev shadow) region is shown in Fig. QJ,. 

Let us find the classic action variation first for trajecto- 
ries with impact parameters b > b* or b < 0, when there 
is no Andreev reflection and the incident roton stays at 
the same branch after the collision. Taking into account 
variation of the roton group velocity along the trajectory, 
Eq. (|T5)) , the variation of the action along the trajectory 
for the incident momentum p > po is 



5S(b)= / (p(y)-p)dy 




1 dy, 



(20) 



F(m) 



E{m) 



r/2 



d9 



o yl — m sin 2 9 

7T/2 

\J\-m sin 2 6d6 



(21) 



ii 



are complete elliptic integrals of the first and the second 
order respectively. In the limits b — > ±oo Eq. (l2Cfl) reduces 
to Eq. (|T5j) . 

In the interval < b < b* trajectory ends at the An- 
dreev reflection point with the coordinate y = —y*. The 
incident roton with momentum p = po + (p — po) > po 
returns after the Andreev reflection to y = — oo at the 
other branch with the same energy but a slightly differ- 
ent momentum p_ = po ~ (p ~ Po) < Po- The variation 
of the action along the whole path is 



SS(b) = / p(y)dy+ / p_(y)dy- / pdy- / p_ dy 




= 2(p-po) 



Here a is an undefined constant, which does not depend 
on b and therefore has no effect on the scattering angle 
(p. Choosing a = one eliminates any discontinuity of 
S(b) at b = and b = b* 2^ Introducing the angle variable 
again one obtains the expression 



SS(b) = 2 J p(y)dy- 



pdy 



= 2(p - po) 



(p-Po) 



(b-b*)F[q> 




b*b 



b 2 + y 2 



l\dy 



-bE 



(23) 



in terms of incomplete elliptic integrals 

de 



F((b, m) 



o VI 



E((f>,m) = / Vl 
Jo 



(24) 



where <fi = arcsin \Jbjb* . 

In Figs.QJi andQJ the action 5S(b) and the scattering 
angle <p(b) — dSS(b)/db are plotted as functions of the 
impact parameter b (in dimensionless variables). The 
angle <p has weak singularities at b — and b = b*, which 
are integrable in the integral for the transverse cross- 
section <7j_ [Eq. (flTj) ]. 

Substituting the transverse cross-section (fTo]) into the 
expression for the parameter D' , which determines the 
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Iordanskii force [see Eq. J3J], one obtains D' = —Kmn n , 
where n n is the normal particle density. This rather sim- 
ple and universal expression tempts to claim its universal 
topological origin, since n in this expression is a topolog- 
ical charge. However, in the Sec. IHII we shall see that the 
expression is not universal, in fact. For quasiparticles in 
a BCS superconductor with energy much exceeding the 
gap an additional small factor appears in this expression. 



III. SCATTERING OF BULK FREE BCS 
QUASIPARTICLES BY A VORTEX IN FERMI 
SUPERFLUIDS 

A. The Bogolyubov-de Gennes equations 

The wave function of quasiparticles in the BCS theory 
is a spinor with two components, 



4>{R) 



u{R) 
v(R) 



(25) 



which are determined from the Bogolyubov-de Gennes 
equations^ 2 -: 

h 2 ,. 



2m 



{R) + A{R)e l6{R ^v{R) = eu(R) 



(26) 



— (V 2 + k 2 F ) v(R) + A{R)e- ie ^u(R) = ev(R) . 
2m 

(27) 

Here k F is the Fermi wave number, and the gap A(_R) 
varies in space. The equations correspond to the Hamil- 
tonian with the density 



^ = £ (|VW|2 - 4H2) -£ (|VU|2 



fc>| 2 ) 

A{R)e w(R) u*v + A(r)e~ je(K Vu.(28) 



If a supcrfluid is at rest the order parameter phase 9 
is a constant and the solution the Bogolyubov-de Gennes 
equations is the plane wave 



where 



*'() 



v e 
I 
\ 



X) 



Ah-R 



(29) 



1 - 



(30) 



/ 



The energy is given by the well-known BCS quasiparti- 
cle spectrum e = ± \/£ 2 + A 2 . Here £ = (h 2 /2m)(k 2 — 
hp) f=a Hv F {k — k F ) is the quasiparticle energy in the nor- 
mal Fermi-liquid and v F = kp/m is the Fermi velocity. 
The two wave numbers 



k± 



k 2 p±2Je 2 -A 2 



(31) 



correspond to the particle-like (+) and the hole-like (-) 
branches of the spectrum. 

The Bogolyubov-de Gennes equations are written for 
the wave function of quasiparticles, and, as in the case of 
the Schrodinger equation, there is the continuity equation 



for the probability \u\ 2 ■ 
some point in space: 



to find the quasiparticle in 



where 



ih 
2m 



d(M 2 + H 2 ) 

dt 



(u*Vu - uVu* 



V-9, 



ih 
2m 



[v*Vv - vVv* 



(32) 



(33) 



is the probability flux. Equation (|32[) is a manifestation 
of the conservation law for the number of quasiparticles. 
But the number of quasiparticles and the probability flux 
g are not the same as the number of particles (charge) 
and the particle current j . The Hamiltonian of Eq. (|28|) 
is not gauge- invariant and there is no conservation law for 
the particle number. The Bogolyubov-de Gennes equa- 
tions lead to the following equation for time variation of 
the particle density \u\ 2 



2. 



0(|t 



M 2 ) 



dt 



= -V- j + 2iA{e- i9 v*u- e l8 vu*), (34) 



where 



\U Vll-UVM (V V!) 

2m y ' 2m y 



vVv* 



(35) 



is the particle current. Equation (IM)) includes sources 
(the last term in the right-hand side) related with possi- 
ble changing of the total particle number. Globally the 
number of particles is of course a conserved quantity. The 
sources in the continuity equation for the particle density 
correspond to conversion of the superfluid part of the liq- 
uid to the normal one and vice versa in inhomogeneous 
states. In order to restore the global conservation law 
one should solve the Bogolyubov-de Gennes equations 
together with the self-consistency equation for the order 
paramour proportional to the gap. This property of the 
Bogolyubov-de Gennes equations is well known in the 
theory of superconductivity 23 . 



B. Superfluid motion in the Bogolyubov-de Gennes 
equations 

Superfluid velocity is determined by the order param- 
eter phase gradient: 



2tt 



-V0, 



(36) 



where k c — h/2m is the circulation quantum for the 
Cooper-pair condensate and m is the particle mass. As- 
suming constant absolute value of the gap A and gradi- 
ent of phase, the solution of the Bogolyubov-de Gennes 
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equations is 



u e 
v e 



i(fe+V8i)-H 
.(fc-V6» 2 )H 



(37) 



Here we introduced separate phases 8\ and 82 for two 
spinor components. Their sum determines the order pa- 
rameter phase 6 = 81 + 62 ■ The spinor ([37)) corresponds 
to the energy (neglecting terms of the second order in 
phase gradients) 

r, \ , hK„ V8 de Vgi - V8 2 



2tt 



dk 



eo(fc) + Hk-v s 



A 2 y gl _ V 2 
fc . (38) 



It looks as if the phase difference 8\ — 82 were of no im- 
portance since it can be removed by redefinition of the 
wave vector fc. Choosing 8\ — 82 one obtains the ex- 
pression for the quasiparticle energy following from the 
Galilean invariance and well known from textbooks on 
superconductivity 2 - 2 -: e = eq + hk ■ v s . But another 
choice is required dealing with the theory of quasipar- 
ticle scattering by a vortex: either 8\ — or 82 = 0. 
This is dictated by cyclic boundary conditions for spinor 
components on the closed path around the vortex (see 
Sees. HHUl and [IIlDj. 

For the choice 8\ = 8 2 Eq. (|55|) yields the following 
expression for the current in the plane- wave state: 



hk 

2m 



N(k)v s 



(39) 



So the superfluid velocity contribution is proportional to 



the charge N(k) 



\uo\ 



2 in the state. 



C. Scattering of free BCS quasiparticles by a 
vortex: simple approach 

The mutual friction force has been calculated for pure 
type II superconductors long ago 24 ' 25 . Since the BSC 
theory describes also the superfluid 3 He and the effect of 
the magnetic field is insignificant for mutual friction in 
type II superconductors these calculations are relevant 
also for singular vortices in the superfluid 3 He. In this 
subsection we use a simple approach: geometric optics 
for low energies €q — A <C A and perturbation theory for 
high energies eo 3> A. A more accurate theory based on 
the partial- wave expansion will be considered in the next 
subsection. 

When the energy of the quasiparticles is close to the 
energy gap of the superconductor (£ <C A), the BCS 
quasiparticle spectrum eo ~ A + Vph 2 (k — kp) 2 /2A is 
identical to the roton spectrum with the roton minimum 
momentum replaced by the Fermi momentum hkp and 
the roton mass fj, replaced by A/vp, where vf = hkp/m 
is the Fermi velocity. So the semiclassical theory for 
rotons can be directly applied to such BCS quasiparti- 
cles, and the transverse cross-section for them is given 



by Eq. (|16p . in which the circulation quantum k is re- 
placed by the circulation quantum k c — h/2m for the 
Cooper-pair condensate and the group velocity for the 
BCS quasiparticles is vq — vf£,/cq- Figure [1] illustrat- 
ing scattering of rotons by the vortex is relevant also for 
low-energy quasiparticles scattered by the vortex. The 
phenomenon of the nearly 180% reflection of quasiparti- 
cles from the area of the Andreev shadow shown in the 
figure is important for description of zero-temperature 
superfluid turbulenc o 26 ' 27 . 

If the quasiparticle energy is much larger than the su- 
perconducting gap, the group velocity vg approaches to 
the Fermi velocity vp and the method of classical tra- 
jectories yields the transverse cross-section K c /vp. This 
result does not look reasonable, because the cross-section 
being small still does not vanish in the limit A — > 0. 
Indeed, the partial-wave calculation o 24 ' 25 yielded that 
in the limit of small A/£ the transverse cross-section 
differed from the semiclassical result of Eq. (TTB1) by 
the factor A 2 /2£ 2 . This also followed from the solu- 
tion of the Bogolyubov-de Gennes equations in the Born 
approximation 2 ?, as shown below. 

Let us consider the perturbation theory with respect to 
the gap A and the superfluid velocity v s = (k c /2tt)"V8. 
For the sake of simplicity the wave vector fc lies in the 
plane normal to the vortex axis. In our case the super- 
fluid velocity is the velocity v v induced by the vortex. In 
the zero-order approximation u ~ cxp(zfc • r) and v = 0. 
In the first-order approximation the second Bogolyubov- 
de Gennes equation (|2~T|) yields 



A exp(— it 



Aexp(-i0) h 2 



£(k)+E(k) [£(*) + E(k)} 2 m 



(fc • V0) > e lh - 



Akr 



(40) 



The first term in curled brackets yields a correction to 
the quasiparticle energy cx A 2 , but does not contributes 
to scattering which is determined by the order-parameter 
phase gradients. So we keep only the second term pro- 
portional to V#. Inserting it to the first Bogolyubov-de 
Gennes equation (|26|) one obtains the following equation 
for the correction u' to the quasiparticle amplitude u ~ 1: 

(V 2 + fc 2 K - (fc- V8)^e lkr . (41) 

This equation is similar to the wave equation for the 
sound wave propagating pass the vortex 4 * 2 ^ and using 
this analogy one easily obtains the expression for the 
transverse cross-section: 



o± = 



A 2 7T _ A 2 K c 

2^k^ ~ Wvf' 



(42) 



The cross-section vanishes at A — » as expected. But 
the question where the geometric optics went wrong still 
remains. The answer is that the cyclic boundary condi- 
tions were violated with the choice 9\ = 8 2 = 8/2. Let 
us move the spinor given by Eq. (|37[) along a closed path 
around the vortex line. After closing the path the phase 
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9 obtains the shift 2-k but the shifts of the phases 9\ and 
02 are equal to n. So the periodic boundary conditions 
for the spinor components u and v are violated. They 
are satisfied only if either 8\ or 82 vanishes. According 
to Eq. (f38|) this modifies the expression for the quasipar- 
ticle energy in the vortex velocity field: 



e = e (k) + (hk ± mva) 



(43) 



Then the value of p in Eq. (|15|) must be replaced by hk ± 
ravQ. Choosing - for quasiparticles and + for quasiholes 
(this is dictated by a physically reasonable condition that 
the cross-section vanishes far from the Fermi surface) one 
obtains the transverse cross-sectior 



,24,25 




(44) 



cr_L 



In the limit eo 3> A this agrees with the expression (j42|) 
obtained from the perturbation theory with respect to A. 
A more rigorous partial-wave analysis of the next subsec- 
tion confirms this result for any ratio A/eo . This provides 
an explanation for shortcoming of the naive geometric- 
optics analysis: It ignored peculiarities of the Aharonov- 
Bohm effect for BCS quasiparticles and used an improper 
definition for the quasiparticle phase shift along the tra- 
jectory. We shall continue the discussion of this issue in 
the end of the next subsection. 



D. Partial-wave analysis of scattering of free BCS 
quasiparticles by a vortex 



The partial-wave analysis in the cylindric coordinates 
r, </>, z uses expansion of the spinor components in eigen- 
functions e ^ of the orbital moment. In the presence of 
a vortex the phase of the order parameter Ae around 
the vortex is 9 = <f> and the partial wave expansion for 
the wave function is 



E 



il<j> 



E 



i{l~l)4> 



(45) 



where ui and vi must satisfy the Bogolyubov-de Gennes 
equations for partial waves: 



equations for partial waves: 

h 2 f d 2 m ^ 1 du t l 2 ui 
2m \ dr 2 r dr r 2 
d 2 vi | Idvi {l-l) 2 vi 
dr 2 



2m 



r dr 



Am = 
-A Ul = 



h 2 k 2 F 

2m 

h 2 k 2 F 

2m 



vi, 



Vi. 



(46) 



In our case the orbital number I is not an ideal quantum 
number since the two components of the spinor corre- 
spond to two different orbital numbers I and I — 1. 

In order to find the scattering phases, we shall look for 
the semiclassical solution of the Bogolyubov-de Gennes 
equations for the scaled amplitudes Ui = uiy/r and Vi = 



viJr: 



d 2 Ui I 2 -1/4 



2m \ dr 2 



Ui + AVf = e 



h 2 k 2 F 

2m 



d 2 Vi (Z - l) 2 - 1/4 1 



2m V dr 2 



Vi)+MJi 



2m ) 1 
(47) 



The semiclassical solution of the Bogolyubov-de Gennes 
equations (|47p for partial waves is 




Jf r k±(r)dr 



(48) 



where e Q = e - (I - l/2)/2r 2 and 



{1-1/2Y 



±2\ 



-A*. 



2r 2 



(49) 



If a quasiparticle with the wave number k + is incident 
on the vortex line, it will be reflected either as a quasi- 
particle with the same number k + (usual reflection) or 
as a quasiparticle belonging to the hole-like branch with 
k- < kp (Andreev reflection). The usual reflection oc- 
curs at the turning point determined by the condition 
fc-l_(rv) =0. In the Andreev reflection point the inner 



1-1/2 
TFT 



± A = 0. The 



radical in Eq. (|4"9"|) vanishes, i.e., e • 
type of the reflection depends on which turning point is 
reached earlier: usual or Andreev reflections take place 
if r r > r a or r r < r a respectively. 

In the following we shall look for the wave function for 
large orbital numbers I, which correspond to large impact 
parameters. Then only usual reflection is possible, and 
one can expand the inner radical in Eq. (|49[) with respect 
to (I - 1/2) /r 2 . Then 



k 2 w k 2 _ (l-l/2) 2 ±(l-l/2)e/Ve^^ 



(I - l/2±e/2Ve 2 - A 2 ) 2 



, (50) 



where k± are values of k±(r) at r — > 00. The total phase 
accumulated after quasiparticle motion from very large r 
to the turning point and back to large r is 



$1 = 2 




k 2 ± d-l/2±e/2Ve^^) 2 dr 



--=2k ± r- 



2 2Ve 2 - A 2 



2 (51) 



Here the phase shift — ir/2 originates from the close vicin- 
ity of the reflection point where the semiclassical ap- 
proach becomes invalid. 2 - In order to find the phase shift 
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from scattering of the quasiparticle (particle branch, the 
upper sign in the expressions above) by the vortex one 
should subtract the phase shift $;o = 2k + r — + 1/2) 
of the /"-partial wave function in the uniform state. This 
follows from the expansion of a plane wave in partial 
waves. Then the scattering phase shift is 



5, = 



2 2Ve 2 - A 2 
7r / e 
~ 4 V 1 ~ Ve 2 - A 2 



signZ.(52) 



The variation of the classical action along the quasiparti- 
cle trajectory is connected with the quantum-mechanical 
scattering phase shift by the relation SS(b) = 2M;, where 
b « I /hp. Thus one obtains that 



AS(±oo) = - 



KkpK c 



1 

VQ 



1 

V F 



(53) 



Inserting it into Eq. ([14]) yields the transverse cross- 
section Eq. (|44|) obtained after correction of the 
geometric-optics expression. In the case of the hole 
branch (the lower sign in the expressions above) one 
should subtract from the phase shift $; in the vortex 
state the phase shift 3>(/_i)o = 2/j_r — tt(\1 — 1| + 1/2) of 
the (Z — Impartial wave function in the uniform state. 

The second term oc 1/vf in the transverse cross-section 
(1441) was considered as anomalous and interpreted in the 
terms of spectral flow-^^ though its original derivation 
from the partial-wave analysis^! did not used this con- 
cept (see further discussion in Sec. lIX|) . The analysis pre- 
sented here demonstrates that it can be explained within 
the framework of the scattering theory taking into ac- 
count peculiarities of the Aharonov-Bohm effect for BCS 
quasiparticles. 



IV. BOUND VORTEX-CORE STATES AND THE 
CURRENT IN THE CORE 

A. Bound Andreev states in a planar SNS junction 

Though the core states were found accurately long 
ago2£ it is useful for the analysis of their role in vortex dy- 
namics to consider a simplified approach to them based 
on geometric optics. Such an approach was suggested by 
Stone—. He also used the model of a normal vortex core 
exploiting its analogy with the ID problem of Andreev 
bound states in the ballistic Superconductor - Normal 
metal - Superconductor (SNS) junction. We also shall 
investigate this analogy. A number of authors addressed 
in the past the question whether and what Josephson 
current is possible through such a junction in full ab- 
sence of the order parameter in the normal layer— ~—. 
They concluded that the Josephson current is possible 
due to phase coherence of the Andreev states, which are 
sensitive to the phase difference on the junction. 



We consider layers normal to the axis y and look for 
the state with the energy 



e = eo + hk ■ v s « eo + • v s 



(54) 



where the wave vector k^{k x , kf, k z ) has a modulus equal 
to the Fermi wave number kp, so that the component k y 
is equal to kf = y 'k F — k 2 . The wave function, which 
satisfies the Bogolyubov-de Gennes equations, is given 
by 



A Am(v s ■ R) j h+imeoy / h 2 k f 



g — im(v s ■ R) /h—imeoy /h kf 

inside the normal layer < y < L, 



e lk °- H (55) 



u _ e ie-/2+im(v s -R)/h \ iko . R+ ^£z3. y 



-i0_/2— im(v s -R)/h 



(56) 



inside the superconductor at y < 0, and 



u \ _ ( u+e iO + /2+im(v s -R)/n \ iko . B ^ m V* k ~ y 
v ) - ° [ v+e -ie+/2-im{v a -R)/h I ' 



inside the superconductor at y > L. Here 



u± 



\ 




(57) 



(58) 



the constants A and B are determined from the boundary 
conditions at the interface y = 0, 



-16-/2 



(59) 



and one can find the value of the constant C satisfying the 
boundary conditions at the interface y — L for discrete 
values of the energy eo satisfying the following Bohr- 
Sommcrfeld condition^: 



2meoL 



h 2 k 



= 2tt 



/ 



1 

S+ 2 



-6L)-2arcsin^|, (60) 



with integer s. At small energy eo <C A this yields the 
spectrum of the Andreev bound states: 



(2mL 



\h 2 ki 



2n 



(61) 



This form of the wave function assumes that only the 
Andreev reflection occurs at the core boundaries, so the 
wave vector is always close to the Fermi surface and its 
normal component k y w kf does not change a sign at 
the reflection. This is a valid assumption in the weak- 
coupling limit when the superconducting gap A is small 
compared to the Fermi energy ep = h 2 k F /2m. In this 
approximation the boundary conditions at the interfaces 
y = and y = L require continuity only of two spinor 



No phase difference across 
the normal layer, the superfluid 
moves outside the normal layer 



Constant superfluid phase Phase difference across 
over the whole bulk the normal layer, no superfluid 
motion outside the normal layer 



FIG. 2. Bound states inside the superconducting gap. Occupied levels below the Fermi level e = are shown by black circles 
while empty levels above the Fermi level are shown by white ones. Center: No phase difference across the normal layer, no 
superfluid velocity inside superconducting areas. Right: There is phase difference 9+ — 6>_ across the normal layer but still no 
superfluid velocity inside superconducting areas. Arrowed dashed lines show shift of level relatively to the gap and the Fermi 
level. Some levels exit from the gap at the top of the gap while some new levels enter the gap at the gap bottom. Left: No 
phase difference across the normal layer, but there is the superfluid velocity v 3 at the bulk of superconductors. This shifts the 
gap with respect to the Fermi level and the number of occupied levels changes. 



components, ignoring the continuity conditions for spinor 
gradients. It is worthwhile of noting that this approxima- 
tion does not provide exact continuity of the quasiparticlc 
current at the interfaces as required by the the conser- 
vation law for quasiparticles. This can be achieved only 
in the next approximation taking into account the pos- 
sibility of usual reflection changing the direction of the 
momentum normal to the layers. However, for our goals 
we need not the quasiparticle current but the genuine 
particle current, and the used approximation is sufficient. 

Now let us find the contribution of bound states to 
the momentum of the liquid. Since in the bound states 
particle and hole components are equal in amplitude, the 



charge N(k) 



in these states vanishes, and 



according to Eq. ([55)) the current normal to the layers 
in any bound state is about Hkf/m. Let us consider the 
case of a wide normal layer L 3> hkp/mA when there is a 
large number of bound states and the sum over them can 
be replaced by an integral. Then the total contribution 
of the bound states to the current is simply a difference 
of the number of states with the wave vectors in opposite 



directions. At T = all states with e = eq + hko ■ v s < 
are filled. This yields that the total momentum in the 
bound states (per unit area of the SNS junction) normal 
to the layers is 



Pbs = - 



h 2 kjv s dfcy 
5e An 2 



(62) 



k<k f 



where fen is the component of the wave vector in the layer 
plane, n is the particle density and 



5e = 7r 



mL 



(63) 



is the distance between discrete levels. The momentum 
corresponds to the current density jb s = Pbs/Lm = —nv s 
inside the normal layer. This shows that quasiparticles 
in core bound states play a role of the effective normal 
component, which exists even at zero temperature^. 



Note that the the phase difference 9 + — 0_ between the 
superconducting banks of the SNS junction has no effect 
on the current in the continuum limit. This is because 
the main effect of the phase difference is a shift of the 
levels in the forbidden gap but not an essential change of 



their number. Indeed, the shift of the levels leads to an 
entry of a new level on one side of the gap and an exit of 
an old level on another side (Fig. [2]). Normally these two 
events are not synchronized and the number of level can 
fluctuate with ±1 level. This fluctuation is not essential 
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in the limit of large number of levels, i.e. in the case of 
not too small superfluid velocities. On the other hand, 
exactly this small variation of the level number leads to 
the Josephson effect in the SNS junctional - — , which is 
beyond the scope of the present work. 

The shift of levels constitutes the phenomenon of spec- 
tral flow, which arises if the phase difference 8 + — 9- 
monotonously increases or decreases and the levels cross 
the forbidden gap. On the basis of analogy between the 
vortex core and the SNS junction the spectral flow was 
believed to exist in the core of the moving vortex either. 
Later we shall see why the aforementioned analogy does 
not go so far. 

One should also consider the contribution of the con- 
tinuum delocalized states with negative energy, which are 
fully occupied in the ground state. For a delocalized state 
(| eo | > A) of a quasiparticle propagating from y = — oo 
to y = co, the spinor in the normal layer < y < L is 
given by the same expression as Eq. (|55[) for the bound 
state, whereas in superconducting layers the states are 
described by spinors 



u e 
v e 



i6 + /2+im(v 3 -R)/h 
-i0 + /2—im(v 3 -R)/H 



iko-R-\-i- 



(64) 



for y > L and 



+r 



u e 
v e 



i0-/2+im(v s -R)/h 
-i0-/2—im(v 3 -R)/h 



^my/el-A 2 y/h?k f 



v e 



i6- /2+im{v s -R)/h 



U e 



— i9- /2-im(v s -R) /h 



-im-s l /eg—A 2 y/h 2 kf 



(65) 



for y < 0. Here t and r are amplitudes of transmis- 



sion and reflection (|i| 



1) which are determined 



from the continuity of spinor components (but not their 
derivatives!)^ at y — and y — L. As well as for bound 
states, the analysis considers only the Andreev reflection 
neglecting probability of usual reflection, which change 
the direction of the wave vector. The amplitudes of the 
spinor components in the normal layer [see Eq. f|55[t] are 



B = tv e~ i8+/2++im( 'V' e °~ A2+eo)L/h2kf . 
The transmission probability is 

f 2 - A 2 

11 e 2 , - A 2 cos 2 [e mL/h 2 k f ~ (6+ - 0_)/2] 



(66) 



(67) 



The transmission probability differs from unity in the 
small energy interval of the order A. So the effect of re- 
flection is not essential for the contribution of delocalized 
states to the supercurrent, which can be found by sum- 
mation of the Eq. (f39| over the whole continuum of free 
bulk states. The whole particle density is accumulated 
in delocalized but not bound states. Neglecting reflection 
for the continuum states, the density and the current in 



the normal and the superconducting areas do not differ 
essentially. So the whole ensemble of delocalized quasi- 
particles is a liquid of nearly constant density n mov- 
ing with the spatially uniform velocity z; s .This demon- 
strates ideal transparency of the ballistic normal layers 
for the supercurrent of delocalized quasiparticles. Note 
that scattering of continuum states by impurities is im- 
possible since all continuum states are fully occupied. 

Summing the momenta in localized and delocalized 
states inside the normal layer, one obtains that the total 
momentum and the current eventually vanish there (with 
accuracy of the small parameter of weak coupling A/e^). 
Keeping in mind the presence of the current nv s in su- 
perconducting layers, this violates the conservation law 
for the particle number, since backflow in our ID geome- 
try is impossible and the current must be constant along 
the direction normal to the layers. A proper conclusion 
from this is that the superfluid transport (but not dif- 
fusive transport with dissipation!) with high superfluid 
velocity in this one-dimensional geometry is impossible. 
But this does not rules out the superfluid transport with 
very low superfluid velocities v s < H/mL when discrete- 
ness of the Andreev bound states and the phase difference 
across the normal layer cannot be ignored. This returns 
us again to the problem of the Josephson effect in the 
SNS junction^i.^. 



B. Bound vortex-core states: ballistic normal core 

Now let us consider bound states in a normal core of 
a vortex. A reliable assumption is that a quasiparticle 
inside the core, where the superconducting order param- 
eter vanishes, moves along an approximately straight tra- 
jectory back and forth reversing its direction of motion 
via Andreev reflection at the boundary of the core. The 
trajectory is chosen to be parallel to the y axis. For tra- 
jectories with impact parameters much less than the core 
radius the bound states are similar to those in the SNS 
junction with the normal-layer width L equal to the core 
diameter 2r c . On the other hand the phase difference 
6 + — 8- = 8 V + 9 S consists from the phase difference in- 
duced by the vortex, 9 V — tt — 2 arcsin(6/r c ) w tt — 26/r c , 
and the phase difference 9 S produced by the superflow 
past the vortex. Here b — M/kt is the impact parame- 
ter and I is the quantum number of the discrete angular 
momentum. Geometry of the process is shown in Fig. [3] 

Eventually the energy of the bound state in the normal 
core for the chiral zero-crossing branch s = — 1 depends 
on the orbital quantum number I and is e(l) = £o(l) + 
hkfV s i cos a, where 



eo(0 



2mr r A 



2mr r 



(68) 



and a is the angle between the trajectory (the axis y) 
and the local superfluid velocity v s i, which is different 
from the superfluid velocity v s far from the vortex in the 
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FIG. 3. (Color online) Bound state in the normal core. The 
vertical solid arrowed line shows the trajectory of the quasi- 
particle and the vertical dashed arrowed line shows the tra- 
jectory of the quasihole after Andreev reflection of the quasi- 
particle at the core boundary. Note that the picture is purely 
schematic, and in fact the analysis was done for the case when 
the impact parameter 6 is much less than the core radius r c 
but still much larger than the interatomic distance 1//cf. 



presence of the backflow (see below). Introducing the 
angular momentum L z = hi of the bound state, there is 
a frequency 



UJQ 



de 
dL z 



2mr? 



2mAr 



(69) 



with which the trajectory slowly rotates around the vor- 
tex axis. The phase difference from the superflow outside 
the core can be presented in the form of a dipole held 
6 S = (1 + h 2 kf /2mr c A)('im/h)r c ■ v$, where r c is the 
vector of the modulus r c directed normally to the cylin- 
dric border of the core and the superfluid velocity vg is 
one more superfluid velocity different in general from the 
asymptotic velocity v s and the local velocity v s i. Intro- 
ducing the isotropic part of the spectrum 



the spectrum becomes 



e = e 00 + hk- (v s - v L ), 



(70) 



(71) 



where the presence of the vortex velocity points out 
that the calculation is done for the coordinate frame con- 
nected with vortex and 



v s = v g + v si 



(72) 



is the effective superfluid velocity taking into account the 
superfluid velocity v s i outside the core and the contribu- 
tion of the phase difference 9 S . 



For the analysis of the vortex mass one needs to know 
the contribution of bound states to the total momen- 
tum of the liquid if the fluid flows past the vortex with 
the superfluid velocity v s . As well as in the case of the 
SNS junction, every bound state has a momentum of 
the magnitude about Hkf directed along the bound-state 
trajectory. Taking into account that the energy interval 
between levels is 6e = Hooq and integrating over all di- 
rections and the wave number component k z along the 
vortex axis one obtains that the total momentum is 



Pbs = - 



1 f dk z f , 2 hk 2 f 

- / — / da cos a — -(v s i-v L 

2 J 2vr J uj q 

-k F 



= [Vsi - v L ). (73) 

UJQ 

The expression was derived for large r c neglecting cor- 
rections of the order h 2 kf /2mAr c in Eq. (|69[) . The 
coefficient before the velocity is the Kopnin vortex 
mass fix — Trhn/ujQ. Using Eq. (|69p for wo in the limit 
of large core r c h 2 kp/mA and assuming that the 
momentum is uniformly distributed over the area nr 2 
of the core this momentum corresponds to the current 
jb s = Pbs/^r 2 m = —2n(v s i — vl) inside the core. 

The momentum in the core depends only on the local 
superfluid velocity v s i outside the core and not on the 
phase difference 8 S . The reason is the same as in the 
case of the SNS junction: the phase difference shifts en- 
ergy levels but does not change their number since any 
crossing of the zero energy by a level is compensated by 
an entry or an exit of an level at the bottom of the forbid- 
den gap. Meanwhile, it is variation of the phase difference 
what leads to the spectral flow. So the spectral flow has 
no effect on the total momentum in the core, and on the 
vortex mass as a result of it (see below). In contrast to 
the SNS case, the spectral flow in the core of the moving 
vortex can only oscillate but not monotonously cross the 
forbidden gap. The oscillation is related with rotation 
of the bound state with angular velocity ujq and depen- 
dence of the level position with respect to the gap on the 
a-dependent phase difference 9 S in Eq. (|6"5)) ^£ 

In the model of the normal core the Andreev reflection 
for all bound states occurs at the core boundary, and 
the energy of bound states are easily determined analyti- 
cally from the semiclassical Bohr-Sommcrfcld condition. 
Meanwhile, the more realistic model of the core with lin- 
ear growth of the order parameter in the core considered 
in App. |A1 shows that though the concept of well defined 
trajectory (geometric optics) works well, one cannot use 
the semiclassical approach for description of motion along 
the trajectory and the Bohr-Sommerfeld condition is in- 
valid. Despite this, the model of the normal core gives a 
qualitatively correct energy spectrum, different from that 
from more accurate theories only by a numerical factor. 
On the other hand, this model allows a simple analytical 
analysis of the backflow effect on the vortex mass, which 
would require less transparent numerical calculations in 
more realistic models. 
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V. VORTEX MASS IN THE BOSE LIQUID: 
BACKFLOW AND COMPRESSIBILITY 
CONTRIBUTIONS 

In an ideal liquid a singular vortex line has no own iner- 
tia and cannot move with respect to the liquid, in which 
the vortex line is immersed (Hclmholtz's theorem). But 
this statement is exact only in the limit of an infinitely 
thin line. Taking into account the finite size of the vortex 
core the vortex line can move with its own velocity v^ 
different from that of the surrounding liquid and there is 
an inertial force proportional to the vortex line accelera- 
tion dvL/dt. 

A naive estimation for the vortex mass is to deduce it 
from the picture of a cylinder moving through a perfect 
fluid assuming that the cylinder has a radius equal to the 
core radius.— Then classic hydrodynamics tells that the 
cylinder induces a dipole velocity field around it (back- 
flow): 



Vh 



2tt 



V0, 



bf 



v bf 



(74) 



The condition of the absence of the radial flow through 
the core boundary in the coordinate frame moving the 
vortex velocity vi requires that the constant velocity Vb f , 
which determines the backflow, is Vf,f = i>l — v s . Here 
v s is the superfluid velocity far from the vortex core. We 
consider the case of T = when n s = n. The kinetic 
energy of the backflow is given by 



(vl — v s ) 2 mnr\ 



dr 2 



(v L -v s )-r 



r>r c 



2 ( V L ~ Vs) 2 

nr c mn . (75) 



So this yields the vortex mass (i v equal to fi core = irr 2 mn, 
which is a mass per unit length of the liquid inside a 
cylinder of the radius equal to the core radius r<r^. Later 
we shall call u core a core mass (in contrast to a more 
general term vortex mass taking into account all possible 
contributions to the mass of the vortex). 

The vortex mass can be determined from calculation 
of the vortex-velocity dependent contribution to the en- 
ergy or the momentum. Naturally the both calculations 
should yield the same mass. Sometimes it is simpler to 
calculate the momentum-. But calculation of the mo- 
mentum of the backflow has a subtlety well known in 
classical hydrodynamics. The direct way to estimate the 
momentum of the potential velocity field in an incom- 
pressible liquid is to integrate the expression for the mo- 
mentum by parts. For the backflow this yields 



P = mn — / 'V9bf(r)dr 



mn— | J 8bfTidS 



ObfndS 



(76) 



where integration is reduced to integrals over the cylin- 
dric surfaces of radius r c and of infinite radius and n 
is a normal to these surfaces. Strictly mathematically 
speaking this yields zero since surface integrals do not 
depend on surface radii for the backflow field. However, 
one should take into account that any finite momentum 
P in an incompressible liquid means that the whole liquid 
moves with the velocity P/mnV inversely proportional 
to the volume V. One should take into account this tiny 
velocity simply by deleting the contribution from the dis- 
tant surface. This yields the momentum called in hydro- 
dynamics Kelvin impulse (see Sect. 119 in the textbook 
by Larnb-^): 



K 



mn 



2tt 



QbfndS = [i core Vbf. 



(77) 



In classical hydrodynamics they justify using the Kelvin 
impulse for an object moving through an incompressible 
liquids by considering the momentum transferred to the 
object when making it to move from rest.— But in quan- 
tum hydrodynamics the justification looks simpler. Lo- 
cal perturbations of the velocity field cannot change the 
phase at infinity. So the boundary condition at infinity 
is not vanishing velocity, but vanishing phase, i.e., the 
potential of the velocity field. On the basis of it the in- 
tegral over the distant surface in the expression for the 
momentum should be ignored. 

The calculation of the vortex mass assumed that a 
moving core is impenetrable for the fluid as a real rigid 
cylinder though the cylinder itself has no mass. In real- 
ity the vortex core is not empty and the superfluid will 
flow through the core, thus producing a reduced back- 
flow field.— So our simple calculation provides only an 
upper bound on the vortex mass related to the core. 
For illustration of this effect let us consider the model 
of a partially filled core with constant particle density 
n(l — A) inside characterized by the parameter A < 1. In- 
side the core the liquid moves with the constant velocity 
Vi n , which corresponds to the phase 9i n = 2ir(vi n ■ r)/ k. 
The continuity of the phases 9i n inside the core and the 
phase 9 out = 27r(t> s ■ r)/n + 9bf outside the core together 
with continuity of the radial flow at the core boundary 
yield: 



v s - v bf , v bf = 



A 



A 



(78) 



In the coordinate frame moving with the vortex velocity 
this gives the momentum 



P L = -nr\mn(\ - \)v in + (S - Ttr 2 c )(v s - v L ) 



+P K = nm(v s -v L )[S- nr 



, 2A 
2- A 



(79) 



where S is the whole area occupied by the liquid. In 
order to see the value of the vortex mass one needs to 
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know the momentum in the arbitrary coordinate frame: 
P = P L + mn[S - nr 2 + (1 - \)-Kr 2 c ]v L 



mn S — nri 



2A 
2 - A 



A 2 



mnnr n 



A 



v L . (80) 



The vortex mass /j, v — [2A/(2 — X)]fi core is a factor before 
the vortex velocity vl- If the density suppression An = 
nX in the core is small the vortex mass is quadratic in 
An. 

Strictly speaking the model of constant density in the 
core is not relevant for singular vortices in Bose super- 
fluids with an obligatory zero of the density on a vortex 
axis. Therefore, in App.[B]we derive the vortex mass for 
the Bose superfluid using a more realistic model with lin- 
ear in r growth of the density in the core. On the other 
hand, the model of constant density in the core is rele- 
vant for continuous vortices in the Fermi liquids, namely 
for estimation of the effect of superfluid density suppres- 
sion on the vortex mass. However, this contribution is 
small compared to the effect of bound states in the core 
(see Sec. ED . 

But in the Bose liquid the most important contribution 
to the vortex mass is connected with finite compressibil- 
ity of the liquid. The cross term in the kinetic energy of 
the velocity field v(r) — vl = v v (r) + v s — vl in the co- 
ordinate frame moving with vortex produces the density 
variation in accordance with the Bernoulli law: 

Sn = -mn—v v (r) ■ (v s - v L ) = — ~v v (r) ■ (v s - v L ), 
aP cj. 

(81) 

where dn/dP = l/mc 2 s is the fluid compressibility, P is 
the pressure, and c s is the sound velocity. The density 
variation leads to the energy contribution 12 - 13 
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-,(82) 



where e is the static vortex energy per unit vortex-line 
length. Like the vortex energy, the vortex mass is de- 
termined a logarithmically divergent integral, which has 
to be cut at some hydrodynamic scale R, e. g., the in- 
tervortex distance. In the Bose superfluid, according to 
the Gross-Pitaevskii theory, the core radius r c ~ k/c s is 
also determined by the sound velocity c s and as a con- 
sequence, the compressibility mass is by the logarithmic 
factor larger than the core mass [i core = irr 2 mn. 



VI. VORTEX MASS IN THE FERMI 
SUPERFLUID 

The two contributions to the vortex mass (from the 
backflow and the liquid compressibility) in the Bose su- 
perfluid in principle are relevant also in the Fermi super- 
fluid. However, the compressibility becomes inessential 



in the weak-coupling limit despite a large logarithmic fac- 
tor. The difference with the Bose superfluid is that while 
in the Bose superfluid the sound velocity goes down (com- 
pressibility goes up) in the weak-interaction limit, in the 
Fermi superfluid the sound velocity remains high being 
always of the order of the Fermi velocity. But the most 
important difference between the Bose and the Fermi su- 
perfluids comes from bound core states, which contribute 
not only to the mutual friction force but also to the vor- 
tex mass 2 -. Earlier we derived the momentum Pi> s in the 
ground state in the presence of the superflow past the 
vortex [Eq. (|73|)] ■ The factor before the vortex velocity 
vl in this expression is the Kopnin mass fix — Trhn/coo. 
However the full vortex mass is not reduced to the Kop- 
nin mass. The current j(, s = P(, s / 'irmr 2 in the bound 
states exists only inside the core and must transform to 
the superfluid current outside the core. The latter cur- 
rent forms the backflow velocity field, which must be de- 
termined from the continuity equation for the total fluid. 
As a result, the Kopnin mass will be renormalized by the 
backflow effect. 

In analogy with the analysis of the backflow for the 
Bose liquid, the local superfluid velocity v s i = v s + Vf,t 
at the core boundary differs from the superfluid velocity 
far from the vortex and the continuity of the current at 
the core boundary is 

jbs = — — ~n(v L - v s - v bf ) = nv b f. (83) 

Note that the current in the continuum of delocalized 
states does not affect this condition because it has no 
discontinuity at the core boundary and contributes the 
same term mnv s i to the two sides of this equation. The 
latter yields 



(J-k 



He 



UK 



-{v L - v s ), 



(84) 



and the total momentum including the backflow momen- 
tum (Kelvin impulse) is 



bs + Trr c mnv b f - 



/'-<:: 



(v L -v s ). (85) 



According to this expression the Kopnin mass [Ik is 
renormalized by the factor 2/Lt core /(/^ + [i cor e) equal to 
4/3 for the value of oj given by Eq. (|6T)|) in the limit 
of large core radius r c 3> H 2 kf /2mA. The most impor- 
tant outcome of this analysis is not this numerical factor, 
which depends on the model of the core anyway, but a 
more adequate insight into the origin of the vortex mass. 
If the Kopnin mass [Ik is much smaller than the core 
mass Hcore, the Kopnin mass is renormalized by the fac- 
tor 2, i.e., the backflow gives the same contribution as 
the bare Kopnin mass. The case of small normal den- 
sity of bound states is realized for a vortex with a con- 
tinuous core in superfluid 3 He when the core radius r c 
essentially exceeds the coherence length £ c = hvp / A and 
Mif ~ Pcoreic/fc- Addressing this case, Volovik— ^ ar- 
rived to an incorrect conclusion that the contribution of 
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the backflow to the vortex mass is negligible compared 
to the bare Kopnin mass. The reason for it was that 
Volovik used the condition of continuity of the superfluid 
component (see Eq. (24.16) in his book 6 ), whereas only 
the total particle number of the liquid but not its super- 
fluid part is conserved in the presence of the Andreev 
reflection. In fact, Volovik estimated the backflow effect 
from weak suppression of the superfluid density inside 
the continuous core considered for the Bose liquid in the 
previous section. He ignored the backflow induced by the 
current in bound states. 



VII. BOLTZMANN EQUATION FOR THE 
CORE-STATES QUASIPARTICLES: THE 
KOPNIN-KRAVTSOV FORCE AND THE 
VORTEX MASS 

If there are impurities in superconductors or collisions 
of bound quasiparticles with free bulk quasiparticles in 
superfluid 3 He, the bound states produce not only the 
vortex mass but also the mutual friction force (Kopnin- 
Kravtsov force). In this case one should use the Boltz- 
mann equation^ Let us write the Boltzmann equation in 
the continuum of semiclassical states bound in the core 
and characterized by the two Hamiltonian-conjugated 
quantities "angle a - moment L z ' : 



dj__d^d£_ de df 
dt da dL 



= d£ 

dL z da dt 



(86) 



The collision term in the right-hand side in the 
relaxation-time approximation is 



dt 



f- f n {e,v n ) 



(87) 



col 



It takes into account elastic collisions with impurities in 
superconductors (then v n is the velocity of the crystal 
lattice) or with bulk free quasiparticles in superfluids. 
Here 



1 



1 



1 



+ 1 



is the distribution function for bound states, which are 
in the equilibrium with the normal component. 

The equilibrium distribution function in the collision 
term has a small anisotropic part if the superfluid part 
moves with respect to the normal part of the liquid. This 
is well known property of the Boltzmann equation in 
superconductors 6 - ' 17 ' 36 . Note that Kopnin— used the dif- 
ferent Boltzmann equation, which follows from that used 
in the paper if the superfluid velocity v s is replaced by the 
normal velocity v n — 0. Though this difference does not 
lead to the difference in the Kopnin-Kravtsov force and 
the Kopnin mass, since they do not depend on the relative 
velocity v s — v n , still it is interesting how superfluid trans- 
port manifests itself in the Boltzmann equation. This can 



be important, e.g., for non-stationary phenomena when 
the distribution function varies in time. 

We expand the distribution functions around the 
isotropic equilibrium distribution function /o(eoo)' 



f{p) = /o(eoo) + fi(e,v n ), 



(89) 



The zero-approximation function /o(eoo) is the equilib- 
rium Fermi distribution function equal to /„ at v n — 
v s = vl. The equation for the first-order correction lin- 
ear in the relative velocities is 



hu k ■ [(v s - v L ) x z\— wo^— 



dfi 



de 



da 



fi-hk- (v n - *s)^r 



(90) 



Its solution is 



/i = ~^ h i k ' - v l) 
w rfc • [(v n - v L ) x z] + k ■ (v n - v L ) 



(91) 



1 + C^ T 2 

The total momentum in the vortex-core bound 
states^ is given by 



— kp J^rnzn 



(92) 



Here L™ ax = L + Hk f r c 9 s /2 and L™ in = -L + 
Hkfr c 9 s /2 are maximal and minimal values of the an- 
gular momentum in the bound state, which differ from 
±L = ±A/ojq because of the phase shift 8 S . This may 
look as if the momentum depends on the phase shift 8 S 
contrary to what was calculated for the ground state in 
Sec. HVBl Indeed, the anisotropic part of the distribu- 
tion function f\ obtained from the Boltzmann equation 
depends on the # s -dependent v s given by Eq. (Tf2|) . This is 
a natural result since the Boltzmann equation takes into 
account only events near the Fermi surface, while entries 
and exits of the bound states to and from the forbidden 
gap at the top and at the bottom of the gap are also 
important, as was demonstrated above. In reality these 
events are accounted for with direction-dependent limits 
L max and L min of the integra i m Eq. One may 

change variable in this integral introducing the modified 
angular momentum L' z — L z — hkfr c 9 s /2 so that 



k p 



bs = 



dk z 

4-7T 



da 



2tt 



h - ^ir hk ■ ( v e - v L ) 
de 



(93) 

The second term in brackets cancel the 6 S dependent 
term in fx, and eventually after using the distribution 
function given by Eq. (l9Tj) only the local superfluid veloc- 
ity v s i appears in the expression for the total momentum 
of core bound states. One can use the modified angu- 
lar momentum L' as a new variable instead of L 7 from 
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the very beginning in the Boltzmann equation (|86l) itself 
with the same result: the phase difference 9 S drops out 
from all expressions and the effective superfluid velocity 
v s reduces to the local superfluid velocity v s i outside the 
core. 

In the limit of zero temperature dfo/de = — 6(e) and 
the momentum in the bound states is 



nhn f U) t[(v„ - v L ) x z }+ v n - v L 

Pbs = 1 v L - v s i H — — ^— 

1 + u>i.T 2 



the mutual-friction force is 



0J 



(94) 



The expression reduces to Eq. (|73[) in the limit t — > oo. 
The part of the momentum linear hi determines the 
Kopnin mass taking into account the effect of collisions. 
The mass becomes a tensor: 



irhnr 



1 



U! T -1 
1 UJ T 



(95) 



Kopnin and Vinokur— related the term in the momen- 
tum transverse to the relative velocity v n — vl with the 
transverse vortex mass. This term, however, does not 
lead to a conservative inertial force, which follows from 
some Hamiltonian. It determines a high-frequency cor- 
rection to the dissipative (longitudinal) mutual-friction 
force, which has its counterpart in the dissipative func- 
tion (see the next section). 

Repeating the process of renormalization of the Kop- 
nin mass by the backflow effect, one obtains the same 
renormalization factor 2fj, core /(fj,K + (J-core) as obtained 
in the previous section without collisions. In this factor 
the Kopnin mass [Ik = nfrn/ujo is the scalar mass in the 
limit t — y oo. 

In the case of frequent collisions (r -C 1/ujq) the veloc- 
ity vl drops out from the expression f|94[) for the momen- 
tum, and the Kopnin mass vanishes. This is a natural 
result since in this limit the effect of reflections from the 
walls of the core is fully suppressed by frequent collisions 
with impurities or quasiparticles. This does not mean the 
total absence of the vortex mass but its absence in our 
approximation, which neglected the effects of the order 
A/ep. It is worthwhile to note that the small uqt does 
not necessarily invalidate the assumption that the mean- 
free path / of quasiparticles is much longer than the core 
radius mentioned in Introduction. Indeed, since r = l/vp 
and wo ~ h/mr 2 the condition luqt <C 1 reduces to the 
condition l/r c <C ep/A. In the weak-coupling limit ep/A 
is very large so even large I jr c can satisfy this condition. 

The contribution of the bound states to the mutual- 
friction force [see Eq. ([1])] is determined by the momen- 
tum transferred from bound states confined in the vortex 
core to normal quasiparticles or impurities via collisions: 



F r = nhn 



UQTjVn -Vl)~ [(Vn - Vl) X z] 

1 + ujIt 2 



(97) 



The force component transverse to v n —v^ is the Kopnin - 
Kravtsov force. In the limit T — >• and luqt —> this 
force compensates the Magnus force, eliminating the to- 
tal transverse force on the vortex. In this limit it has a 
universal r-independent value, hinting to its topological 
origin. For better understanding its meaning let us derive 
it not from integration of the collision term but from in- 
tegration of the left-hand side of the Boltzmann equation 
(|8"6")) . which corresponds to that of the Liouville equation. 
It contains the divergence of the flow of quasiparticles in 
the phase space {a, L z }. Since we are restricted with the 
case of r — » the solution of the Boltzmann equation is 
the equilibrium distribution function /„ [see Eq. (|88[) ] 
for the quasiparticle gas at rest in the coordinate frame 
connected with the normal component. Integrating by 
parts the Kopnin-Kravtsov force is 



c ~ 2 

de df de df 
da dL z dL z da 

,„,„„ LT a * 
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x k]u} f{L z ) 



(98) 



Here we took into account that de/dL z does not depend 
on a and de/ da does not depend on L z . The first term in 
this expression is the momentum flux through the bound- 
ary of the core and the second term is the momentum 
transfer from the external force driving the vortex at 
the process of the bound state rotation with the angular 
velocity wo = de/dL z . While the isotropic part of the 
distribution function contributes to the first term, only 
the anisotropic part provides the second term. Inserting 
/ = /„ one obtains the Kopnin-Kravtsov force in the 
limit ujqt — » 0. So the origin of the Kopnin-Kravtsov 
force looks clear and does not require a reference to the 
artificial concept of spectral flow. This does not rule out 
its topological nature since the force is proportional to 
the circulation quantum, which is a topological charge. 



VIII. EFFECT OF VORTEX MASS ON VORTEX 
DYNAMICS 



dk z 
2tt 



f da J 




df 




2tt 


dt 
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Substituting fi from Eq. (|9"Tj) the core contribution to 



Taking into account all forces discussed above the gen- 
eral equation describing free motion of the vortex is 



dv L 
dt 



-mn M n [z x v L \ = -7»£-Mjl 



Z X 



dv L 
dt 



(99) 
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Here um is the effective density, which determines the ef- 
fective Magnus force. At zero temperature it varies from 
njvf = n for superclean superconductors down to tim = 
for moderately dirty superconductors. The right-hand 
side of the equation contains two dissipative forces. The 
second of them is connected with the transverse vortex 
mass originated from core bound states^ It determines 
a high-frequency correction to the dissipative (longitudi- 
nal) mutual-friction force, which does not appear in the 
Hamiltonian but has its counterpart in the dissipative 
function. In order to see it let us derive the time varia- 
tion of the kinetic energy of the vortex: 



dE a 

lt = dt[ tlv Y ]= - FD - 



Here the dissipative function is 



dt 



(100) 



(101) 



The contribution of the transverse mass to the dissipative 
function is not positively defined. Therefore, the equa- 
tion of motion as given by Eq. makes sense only 
if the transverse-mass contribution is small compared to 
that of the usual friction force cx 7. 

Without dissipation Eq. describes rotation of the 
vortex with the angular velocity equal to lo v — winn/ fj, v ~ 
«/r^ for the vortex mass (i v ~ mnr 2 . The hydrodynam- 
ical approach is valid only if the vortex moves around a 
circumference of the radius ro, which exceeds the core 
radius r c . Then the linear velocity u} v tq must exceeds 
the value of the critical velocity ~ n/r c . Hardly this ro- 
tation is of practical importance. The frequency w v also 
characterizes the frequency of an a.c. process at which 
the vortex-mass effect can compete with the transverse 
Magnus force. In the Bose liquid the frequency u> v is on 
the order c 2 s /n. A phonon with such a high frequency 
has a wavelength comparable with the core radius. In 
the Fermi liquid the frequency uj v is of the same order 
as the frequency Uq, which determines the distance frwo 
between core energy levels. Though the latter is small in 
comparison with the gap A in the weak-coupling limit, 
the frequency itself is rather high. The vortex mass could 
be of more interest in the cases when the Magnus force 
is suppressed by the Kopnin-Kravtsov force. However, 
impurities, which are the source of the Kopnin-Kravtsov 
force, decrease also the vortex mass. In all, it is not 
simple to reveal the vortex mass in superfluids and su- 
perconductors, though some experimental evidence of the 
vortex mass in superconducting thin films has been re- 
cently reported^ 



IX. DISCUSSION AND CONCLUSIONS 

As was already mentioned, Volovik 16 suggested that 
the transverse force from core states (the Kopnin- 
Kravtsov force) and the part of the transverse force from 



the free quasiparticles in the bulk of the Fermi superfluid 
[the term ~ 1 /vp in the transverse cross-section given by 
Eq. flUJ) ] originate from the spectral flow in the core (see 
also Ch. 25 in his book 6 ). He argued that the process of 
vortex motion in space is accompanied by a steady shift 
of core bound-state levels from negative-energy contin- 
uum to the positive-energy continuum. Any crossing of 
the gap by a bound state leads to transfer of the mo- 
mentum, which leads to the transverse Kopnin-Kravtsov 
force. 

Deriving the spectral flow Volovik considered the an- 
gular momentum L' z — z ■ \{r — (vl — v n )t) x p] around 
the axis, which moves together with the thermal bath 
(normal component). Here r is the position vector with 
the origin on the symmetry axis of the vortex moving 
together with the vortex. Volovik's angular momentum 
varies in time: 



dL' 
~~dt 



x PI- 



(102) 



Since the energy of the bound state is proportional to the 
angular momentum, Volovik concluded that the energy 
levels cross the forbidden gap and the the zero energy 
level with the rate proportional to vl — v n . The prob- 
lem with this argument is that the position of the bound 
state energy with respect to the gap depends on the an- 
gular momentum about the symmetry axis of the vortex 
in the coordinate frame moving together with the vortex. 
Then the angular momentum is conserved and provides 
a good quantum number, which determines the energy of 
the bound state. In any other coordinate frame with the 
reference axis, which does not coincides with the vortex 
axis, the angular momentum is not conserved and is not 
a proper parameter for the bound state. Moreover, de- 
riving Eq. (|102l) , Volovik assumed that the momentum p 
of the bound state does not varies in time. Meanwhile, in 
a genuine bound state the momentum p rotates with the 
angular velocity ujq and vanishes in average. As a result, 
dL' z /dt vanishes also, and the angular momentum deter- 
mined with reference to any axis does not differ from the 
angular momentum around the vortex symmetry axis. 
This is a direct consequence of the theorem of mechan- 
ics, which tells that for a system with vanishing velocity 
of the center of mass the angular momentum does not 
depend on the choice of the reference axis. So the vortex 
motion with respect to the thermal bath does not lead 
to the spectral flow. 

Volovik stressed that his derivation was only for con- 
tinuum limit LJ{)T — > when levels are strongly broad- 
ened. However, in this limit the very concept of the 
spectral flow becomes ambiguous. This is illustrated in 
Fig. |U which shows the effect of the level shift on the 
density of states no(e) for various luqt. Without colli- 
sions (loqt — > 00) the density of states is a chain of very 
narrow peaks and a shift of the levels with respect to the 
forbidden gap is a clear effect (Fig. @Ji). For very small 
but still finite ojqt the effect of level shift on the density 
of states is much weaker but still noticeable (Fig. HJd) . In 
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ing banks monotonously varies in time, as, e.g., in the 
a.c. Josephson effect. But the superfluid phase differ- 
ence across the core of the moving vortex does not vary 
in time. Therefore, the spectral flow exists in the former 
case, but is totally absent in the latter. 

The absence of spectral flow in the core automatically 
rules out the spectral flow in the continuum of delocalized 
states also suggested by Volovik— . Indeed, in stationary 
processes these spectral flows should be equal. Otherwise 
there were accumulation or depletion of states at the bor- 
ders between localized and delocalized states. It is shown 
in this paper that the whole transverse force on the vor- 
tex from delocalized states in Fermi superfluids can be 
explained by peculiarities of the Aharonov-Bohm effect 
for BCS quasiparticles without referring to the spectral 
flow. 

The paper clarifies also the role of the backflow on the 
vortex mass. The backflow is an ubiquitous phenomenon, 
which arises from mismatching of currents inside and out- 
side the vortex core (either due to suppression of the fluid 
density in the Bose liquid, or due to to currents through 
core bound states in the Fermi liquid. Its existence fol- 
lows from the conservation law for the particle number 
(charge) . In the Fermi liquid the backflow leads to renor- 
malization of the Kopnin vortex mass by a numerical fac- 
tor both for singular and continuous vortices. 



FIG. 4. Effect of shift of energy levels on the density of 
states at various uiqt. The density of states is shown by solid 
lines before the shift and by dashed lines after the shift, a) 
ujqt — > oo. The density of states is a chain of sharped peaks, 
b) ujqt <C 1. Very broad peaks strongly overlap and cause 
only weak oscillations of the density of states, which are still 
noticeable in principle, c) ujqt — 0. The plot of the density of 
states is totally flat and its shift does not lead to any physical 
consequence. 
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the extreme case ujqt — when oscillations of the den- 
sity of states are totally undetectable the level shift does 
not lead to any effect and cannot influence any physical 
process. Without taking into account whatever tiny os- 
cillations of the density of states it is impossible even to 
define it. 

Altogether this puts in question not the Kopnin- 
Kravtsov force itself but the connection of the force with 
the spectral flow. So the claim that the experiment on 
mutual friction force confirms the spectral flow^S is not 
justified. It is the Kopnin-Kravtsov force, which was re- 
vealed in the experiment, but not the spectral flow. 

Arguing for the spectral flow in vortex dynamics 
they frequently draw an analogy with the Andrecv- 
reflection bound states in the Superconductor-Normal 
metal-Superconductor (SNS) layered system. Though 
this analogy is useful indeed for the bound states in the 
model of normal coreii it fails with respect to the role of 
spectral flow. In the SNS system the spectral flow really 
exists if the phase difference between two superconduct- 



Appendix A: Bound states in a core with linear 
growth of superfluid density 

A reliable assumption is that a quasiparticle inside the 
core moves along an approximately straight trajectory 
back and forth changing its direction of motion via An- 
dreev reflection. The trajectory is chosen as the y axis. 
One can refer to the Bogolyubov-de Gennes equations in 
the ID case eliminating fast oscillations of the wave func- 
tion by the transformation u — ue lkFV , v — ve lkFV . For 
the sake of simplicity we assume that there no compo- 
nent of the wave vector parallel to the z axis. Neglecting 
second derivatives of u and v the Bogolyubov-de Gennes 
equations are 

_ lhvF ^iyl + ^ e *W«(i (V ) = eu(b,y), 
ay r c 

i% VF ^hA + ^ e i9 ^u(b, y) = ev(b, y). (Al) 
ay r c 

Here r = \Jb 2 + y 2 , and the linear dependence of the 
gap A(r) = Ar/r c on the distance r is assumed. In the 
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absence of superfiuid motion through the core the phase 
9 coincides with the azimuthal angle <f> — arctan(y/6), 
and the Bogolyubov-de Gennes equations become 

du(b,y) A(b + iy)_, 



ihvi 



ihvF 



+ 



dy r, 
dv(b,y) A(b-iy)_ 



v(b,y) = eu(b,y), 



u(b,y)=ev(b,y).(A2) 

ay r c 

The normalized solution of the Bogolyubov-de Gennes 
equations is 



u = —v = 



with the energy equal to 



_ bA_ 
£o — — —UqL z 



(A3) 



(A4) 



Here the length f = V r c6= is the geometric average of the 
core radius r c and the coherence length £ c = hvp/ A with 
all three lengths being of the same order of magnitude, 
L z = hkpb is the angular momentum of the bound state, 
and the frequency 



U 



hkpr c 



(A5) 



gives the angular velocity of slow trajectory rotation 
around the vortex axis, in accordance with the canon- 
ical relation equating the rotation velocity to dW/dL z . 

The energy spectrum given by Eq. (|A4[) insignificantly 
differs from the spectrum obtained in the original paper 30 
and in the book by de Gennes— more accurately using 
the partial-wave analysis and a more realistic variation 
of the gap A(r) in the space. This agreement confirms a 
simple picture of the bound states assuming well defined 
trajectories of quasiparticle motion. However, it is nec- 
essary to stress that though the trajectory is well defined 
in the sense that the impact parameter is well defined, 
the motion along trajectory cannot be described semi- 
classically. In particular, our solution shows that there 
are no well defined Andrccv-rcflcction points. Using the 
semiclassical approach and the Bohr-Sommerfeld condi- 
tion for calculation of energy levels one obtains a totally 
wrong spectrum, which is not linear in the angular mo- 
mentum. So the semiclassical theory of motion along the 
trajectory of the bound state is valid only for the model of 
the totally normal core but not for more realistic models 
with non-zero order parameter in the core. 



Appendix B: Vortex mass of a core with linear 
growth of superfiuid density in the Bose superfiuid 

One may consider a more realistic estimation of the 
vortex mass using the Gross-Pitaevskii theory. Accord- 
ing to the Gross-Pitaevskii theory the density grows lin- 
early with the distance r from the vortex axis. Extrap- 
olating this dependence up to the core radius r c and ap- 
proximating the density outside the core by the constant 



value n, the continuous density distribution is 

n(r) = ( n - & l " < T ° . (Bl) 
w [ n at r > r c 

The liquid mass inside the core, 

mn f 2 

m core = 2TT / rrdr = -n core , (B2) 

r c J 3 
o 

is by the factor 2/3 less than the core mass \x core esti- 
mated under the assumption that the liquid density is 
not suppressed inside the core. If the superfiuid moves 
past the vortex the continuity equation in the coordinate 
frame related to the vortex is 

V[n(r)V9] = n{r)V 2 9 + Vn{r) ■ V9 = 0, (B3) 

where the phase 9 determines the velocity field: v s (r) — 
v l = (/s/27r)V0(r). From symmetry all fields are dipole 
fields, and the phase in the cylindric coordinate system 
is 6(r) = 0(r) cos <fi, where <j> is the azimuthal angle with 
respect to the velocity v s — v^. The one-dimensional 
function 9(r) is determined from the equation: 



df_ 

dr 2 



1 



1 dn(r) 
r n(r) dr 



— — % = Q. 

dr r 2 



(B4) 



Inside the core Eq. (|B4[) yields that 9 cx r a with the ex- 
ponent a = (VE— 1)/2 < 1. This means that the velocity 
(but not the current!) has a weak integrable singularity 
at r = 0. The continuity of the azimuthal component 
of the superfiuid velocity at the core boundary r = r c is 
satisfied by the following phase distribution outside and 
inside the core (the azimuthal angle dependence is omit- 
ted): 



(Vs-Vbf)—^. 



(B5) 



Here v B is the superfiuid velocity far from the vortex in 
the coordinate frame moving with the vortex and Vbf 
is the amplitude of the backflow velocity field given by 
Eq. (|74|) . Continuity of the radial velocity gives the con- 
dition: 

d9in , . d9 ut / \ 

a[v s - v b f) = ^ = v bf + v s . (B6) 



dr 

This yields the relation 



v b f = 



1 



, ■ (B7) 
1 + a 

The total momentum includes the momentum Pi n inside 
the core, the momentum of transport superfiuid veloc- 
ity v s outside the core, and the Kelvin impulse of the 
backflow velocity field outside the core [Eq. ([77)) ]: 

Pl = ~2~^ J n ( r ) Vd ™dr + mn(S - nr 2 c )(v s - v L ) 



r<r c 

+Pk = mn 



S — nr 2 



1-a 



2+ a 1+a 



(v s - v L ), 

(B8) 
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where the superfluid velocity v s was replaced by the rel- 
ative velocity v s — v^. The last step is to transform 
the momentum Pt, in the coordinate frame moving with 
the vortex to the momentum in the arbitrary coordinate 



frame: 



P = Pl + [mn(S - nr%) + m core ]v L 



1 - a 



2+ a 1+a 



v s + VvV L - (B9) 



Here fi v is the vortex mass. Taking into account the value 
a = (\/5 — l)/2 the vortex mass is 



pb v = Trr c mn 



irr r mn 



2^5-4:1 



1-a 1 
1+a ~~ 3 

= 0.139^ cor 
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